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DIGITAL VIDEO: MULTIPLE SUB-NYQUIST CODING 

J.H. Stott, M.A. 

GJ. Phillips, Ph.D., M.A., B.Sc, M.I.E.E. 



1. Introduction 

The use of digital transmission to convey video signals 
has been the subject of many studies, and it is likely that it 
will eventually come into service. However, analogue 
transmission will not be displaced overnight and the two 
methods will co-exist for some time. During this period a 
practical link between two points may be partly analogue 
and partly digital, with perhaps more than one section using 
each type of transmission. 

Thus in designing a system in which an analogue 
signal is translated into digital form and subsequently de- 
coded to analogue form again we must consider the degree 
of impairment produced when the process is repeated. 
For example, a criterion which has often been used is to 
require that acceptable signal quality is maintained after 
four successive coding and decoding operations have been 
performed. 

A novel technique for coding composite PAL video 
signals has been described; 2,3,4 it uses a sub-Nyquist sam- 
pling frequency corresponding to twice the PAL colour 
subcarrier frequency.* The method requires the use of 
comb filters which introduce some degradation of picture 
quality, principally a slight loss of high-frequency diagonal 
resolution. The overall degradation is small and, when the 
video samples are described using 8-bit pulse-code modula- 
tion (p. cm.), the resulting picture quality is still acceptable 
for broadcasting purposes.** The method is attractive 
since the reduction of sampling frequency compared with 
previous systems (often using 3/ sc sampling) gives a corres- 
ponding reduction in bit rate, offering the prospect of 
cheaper transmission. 

We must now consider what happens when this coding 
process is repeated. If the effects of the comb filtering 
stages were cumulative, then the picture quality offered by 
the system after four successive coding and decoding 
operations would clearly be unacceptable for broadcasting 
purposes. The purpose of the work to be described was to 
investigate the problem both in theory and in practice. 

We hope to demonstrate that repeated sub-Nyquist 
sampling, of the type under consideration,*** does not in 
fact introduce any fundamental degradation beyond that 
introduced by the first sampling operation. The coding 
system does, of course, suffer the cumulative effects 
(principally quantising noise) of repeated analogue-to- 



* For System I PAL, colour subcarrier frequency = / sc = 4-43361875 

MHz. 
** Subjective tests described in Ref. 3 give a grading of 1-5 on the 

EBU six-point impairment scale. 
*** Much of the analysis is of general application: under certain 
conditions it applies to other types of sub-Nyquist sampling and 
also to super-Nyquist sampling. 



digital-to-analogue conversion, in common with any other 
digital coding system. Experimental work will also be des- 
cribed, intended to confirm the theoretical results. 



2. The sub-Nyquist sampling technique for PAL 
signals 

The sub-Nyquist sampling technique applied to PAL 
colour television signals is fully described elsewhere, ' ' 
and only a brief description will be given here. 

The working of the system relies on the fact that 
the spectrum of the video signal has gaps wherein little or 
no information is normally present. If we consider a mono- 
chrome signal initially, most of the energy of this video 
signal will be clustered around harmonics of the line scan- 
ning frequency. If it is sampled at a rate/ less than twice 
the highest frequency / m present, (i.e. sub-Nyquist), alias 
components will be generated which will overlap the higher 
frequency part of the signal, (Fig. 1). With an arbitrary 
choice of signal it would not be possible to distinguish the 
wanted and alias components but with a video signal we can 
take advantage of the bunched nature of its spectrum. 

Suppose that a sampling rate f is chosen which is 
equal to (n + Vz) times line scanning frequency, where n is 
an integer. Line frequency harmonics in the input signal 
will now cause the generation after sampling of alias com- 
ponents half-way between line frequency harmonics. Thus 
with most of the input signal clustered around line fre- 
quency harmonics, most of the alias components will be 




f +f 

'o T 'm 



frequency 



Fig. 1 - Spectra of video signal before and after sampling 

(a) Spectrum of video signal before sampling 

(b) Spectrum of sampled video signal for/ >2/ m 

(c) Spectrum of sampled video signal for/ <2/ m 
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Fig. 2 - Overall response of comb filter 

between line frequency harmonics. These may be removed 
by a comb filter having nulls at (n + V 2 ) times line scanning 
frequency over the range of overlap (Fig. 2). We have 
already assumed that little energy is present at these inter- 
mediate frequencies in the original signal, so little infor- 
mation is lost. If a similar comb filter is also used to pre- 
filter the signal before sampling, then those parts of the 
input signal which would cause aliasing at line harmonics 
at the output (which the output filter cannot remove) will 
be removed at the price of a further slight loss of input 
information. 

The degradation caused by the system is a slight loss 
of high-frequency diagonal resolution. 

It may be shown that PAL colour video signals may 
also be handled by this technique ' ' provided that: 

(i) the sampling frequency is exactly 2/ sc (which is 

very close to (h + !4) times line scanning frequency), 

and (ii) sampling is performed in the correct phase relative 
to the colour subcarrier. 

Additionally, certain constraints are placed on the form of 
comb filter employed. 

In respect of luminance information the behaviour of 
the system is as outlined above for monochrome signals. 
It may be shown that with the above provisions the PAL 
colour subcarrier signal is also conveyed without serious 
distortion except for the introduction of some loss of reso- 
lution of horizontal colour boundaries. The total degra- 
dation from a single coding and decoding operation has been 
assessed subjectively, using a system with 8-bit linear p.c.m. 
coding, as grade 1-5 on the EBU 6-point impairment 
scale; this would correspond to roughly grade 4-5 on the 
new CCIR 5-point impairment scale (CCIR Recommenda- 



tion 500, Geneva 1974). This performance is probably 
acceptable for a broadcasting signal transmission circuit. 

The block diagram of the basic system is shown in 
Fig. 3; for clarity the means of generating correctly-phased 
sampling pulses is omitted. 



3. Multiple sub-Nyquist coding 

3.1. General considerations 

The introduction of digital transmission of video 
signals is likely to take place on a piecemeal basis. As 
mentioned in the Introduction, any form of digital coding 
used must therefore be capable of giving acceptable quality 
when repeated several times in succession, a useful criterion 
being the quality obtained after four successive coding 
operations. 

If the sub-Nyquist system described in the previous 
Section is to be considered for use in digital transmission 
the effect of applying it repeatedly must be known. In 
particular, the question arises whether the loss in diagonal 
resolution, etc., will build up with successive coding opera- 
tions. This problem is considered in Sections 3.2 and 3.3, 
and, by an alternative approach, in Section 3.4. It will be 
shown that the degradations caused by comb filtering and 
sub-Nyquist sampling are not in theory cumulative. 

A further consideration is the effect of interposing a 
practical analogue link and practical filters between success- 
ive coding operations. This will be discussed in Section 
3.5. 

Finally, Section 3.6 will deal with the effects of 
repeated quantisation. 

3.2. Analysis of repeated sampling 

In this Section we shall consider the effects of 
repeated sampling of the signal. The process of digital 
coding does, of course, involve both sampling and quantis- 
ing, but the latter will not be considered for the present. 

Fig. 4 shows two systems which will be analysed. 
Fig. 4(a) represents the basic sub-Nyquist system consisting 
of a pre-coding comb filter, sampler and post-coding comb 
filter. The comb filter responses are X 1 {f) and X 2 (f) 
respectively. Fig. 4{b) shows a system where sampling has 
been repeated, but with some filter of arbitrary frequency 
response G(f) placed between the samplers. For example, 
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Fig. 3 - Block diagram of sub-Nyquist coding system. The a.d.c. samples at 2f sc 
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Fig. 4- The systems for analysis 

(a) Single-sampling system 



in a simple repetition of the basic system, Gif) would 
consist of a post-coding comb filter and pre-coding comb 
filter in cascade. By comparing the two systems we may 
determine the extra degradations, if any, introduced by the 
second sampling process. 

Certain simplifying assumptions are made: 

1 lT Q is the sampling fre- 



(ii) 



Gif) = 0,\f\>f o where f 
quency. 



S{f) = 0, |/l > f Q where S(f> is the spectrum of the 
input signal s{t), given by the Fourier transform 
relations: 



Sift 



I 



s(t)e~ iwt dt 



(w = 2irf, this 
notation is used 
throughout the 
report). 



I 



sif) = S{f)J ut df 



(iii) The samplers are assumed to have normalised gains 
and have their outputs low-pass filtered so that, for 
example, S 2 (/) = for \fl >f Q . This is valid since the 
samplers are always followed by filters of response 
Gif) or X 2 (f), both of which have this low-pass 
property. The advantage is that high-order alias com- 
ponents may be neglected at an earlier stage in the 
analysis. 

(iv) G{f) is a real filter so that its impulse response g{t) is 
real. 



{b) Repeated-sampling system 
Single-sampling system, Fig. 4(a) 

The filter responses are known, so that 
S 1 if) = X l if) S{f) 



and 



s n (n=x,(f)s,(f) 



d) 

(2) 



The sampler samples s x (t) at a rate/ ; according to 
classical sampling theory we may write 



S,(fl = 



2> 



(f-«/o> 



which, with the assumptions made in (iii) above, simplifies 
to 

5 a V)=S 1 V)+S 1 (r + / )+S 1 (f-/ ) (3) 

This is the result illustrated in Fig. 1. 

Repeated-sampling system, Fig. Mb) 

We are interested in the additional effect of the 
second sampling operation, that is to say the difference 
between signals s (t) ands s U) as shown in Fig. 4. It will be 
useful to define an 'effective insertion transfer function', 
Y(f), of the arbitrary filter, of response Gif), followed by 
the second sampler. 



Let 



i.e. 



Y[f) 



Y\f) 



'S {f) X 2 (f)S 2 (f) 

_ S A V) 
S.V) 



(4) 
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At this stage the mathematics becomes somewhat 
complicated and, for clarity, further working has been 
relegated to the Appendix, Section 7. Section 7.1 derives 
the equation describing the effect of sampling at instants 
t = t + nT , 7.2 applies this to the repeated-sampling 
system and 7.3 derives an expression for the function Ylf) 
we have defined. The result, for positive / in the range 
</</ , is Equation (A6) of Section 7.3: 



Y{f) = G{f)+e-1 w ° T G{f-f ) 



(5) 



For negative /the second term of the R.H.S. must be 
changed to take account of the fact that G(f — f ) is then 
zero (see Equation (A7) of 7.3). An alternative approach 
is to use the result of 7.3 that Ylf) behaves like a real filter 
and therefore obeys the relation 

Ylf) = Y*(-f) 

To summarise, if we know the response Gtf) of the 
arbitrary filter and the delay t between the operations of 
the two samplers we can specify the effective insertion 
transfer function Ylf) which describes how the overall 
response of the system of Fig. 4(b) differs from that of the 
system of Fig. 4(a). 

3.3. Requirements for no added distortion 

We are particularly interested in the condition where 
the insertion of the combination of G(f) and the second 
sampler causes no added distortion, so that s s (t) is identical 
to s Q (t) except for some delay T. Translated into the fre- 
quency domain this implies 

Y(f) = e~i" T 



(6) 



Substituting this condition into Equation (5) we 
obtain the general solution for no added distortion: 

r' uT = G(f) + e~i"° T G(f - f ) for </</„ (7) 

(The reader may obtain a corresponding equation for 
negative / using Equations (6) and (A7) — the rest of this 
section will consider only positive/arid rely on the assump- 
tion that g(t) is real). 

Having obtained this general result it will be useful 
to consider the more special case where G(f) has uniform 
group delay T. In other words 



G(f)=A(f)e->" T 



(8) 



where A(f) is real, and denotes the amplitude-frequency 
response of the filter. 

Substituting (8) into (5) we obtain 

Y(f)= A(f)e-i" T + e -/"or A(f- f )e~i ( " ~ "° )T 

i.e. Y(f) = e->" T \A{n+e->"° W - T) A(f-f )} (9) 

We note that in the practical system which we are modelling 
the second sampler arranges its phase of sampling using the 



colour subcarrier as a reference so that it samples at the 
peaks and troughs thereof.* t is thereby chosen to 
satisfy 



T = T +nT o 



where T c denotes the phase delay introduced by the net- 
work Of) at/ = / /2, the colour subcarrier frequency. If 
the network has uniform group delay T then obviously 
T c = T. Thus t = T + nT Q so that (9) simplifies to 



-/«rj 



Y{f)=e->" T \A(f)+A(f-f ) 



(10) 



and we may see at once that Ylf) is a pure delay T provided 
that 



A\f\+A[f-f Q ) = 1 



(11) 



This is the simple condition for no additional distortion 
when Glf) has uniform group delay. 

Two important conclusions may now be drawn: 

(i) It may be shown that the system comprising a post- 
coding and pre-coding comb filter in cascade has a frequency 
response which satisfies (11) combined with a uniform 
group delay.** Therefore, no extra degradation is intro- 
duced when the sub-Nyquist sampling arrangement of Fig. 3 
is repeated by joining another in cascade. It follows that an 
indefinite number of such systems may be cascaded with- 
out further degradation. 

(ii) It is not necessary to use a comb filter pair between 
sub-Nyquist sampling operations. Any filter with uniform 
group delay whose frequency response A If) satisfies Equa- 
tion (11) may be used. An endless number of such filters 
may be envisaged, but a set of particular interest comprises 
those with smooth low-pass responses of which an example 
is shown in Fig. 5. The response is antisymmetric about 
/=/ /2(=/ sc ), at which point it has the value 1 / 2 . 

This suggests a method of cascading digital and 
analogue links whereby an intermediate analogue link would 
carry a non-standard video signal. Fig. Q(a) shows the 
possible arrangement. The performance would be identical 
to that given by the conventional sub-Nyquist arrangement 
(Fig. 6(b)) but replaces two comb filters by one antisym- 
metric low-pass filter. Assuming 5-5 MHz bandwidth for 
analogue equipment this would require A(f) of the filter to 
bezerofor/>5-5MHz. This implies A(f) = 1 for/from 
to the 'mirror' frequency of 5-5 MHz, which is 8-87 — 5-5 
= 3-37 MHz. Between 3-37 MHz and 5-5 MHz the response 
is a smooth roll-off antisymmetric about 4-43 MHz, at 
which frequency Alf) = 1 /a. The signal handled by the 
analogue link would be non-standard in that the colour 



* This is the same as is performed by the first sampler which handles 

a signal whose subcarrier phase has been made uniform by the 

pre-coding comb filter (see Ref. 3). 

** The analysis for a pair of simple comb filters (combing over the 

whole band) is given in the Appendix, Section 7.4. The analysis 

for the more complex filters actually used is more protracted. 
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Fig. 6 - Cascaded digital and analogue links 
[a) Using a low-pass filter in the intermediate stage [b) Using comb filters in the intermediate stage 



information would be of uniform phase: this is the same 
as the output of the original pre-coding comb filter. It is 
interesting to note that this technique could be used to 
transmit a PAL video signal in a bandwidth of less than 
5-5 MHz (see Section 3.5). 

The possible advantage of such a system is that the 
low-pass filter of the type specified might be cheaper to 
build than a pair of comb filters; against this must be 
weighed the disadvantage that the signal carried on the 
analogue link is non-standard, possibly presenting difficul- 
ties of monitoring. 

A slight modification may be proposed whereby low- 
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pass filters are placed before and after the intermediate 
analogue link. In this case the combined response of the 
two filters should have the form prescribed for A if) above. 
This system would reduce slightly the effects of noise added 
by the analogue link. 



In conclusion, we have shown that it is possible in 
principle to repeat sub-Nyquist sampling without any added 
degradation, provided certain conditions are met. More- 
over, the equations derived may be used to calculate the 
degradations introduced when the performance of a practi- 
cal system departs from the ideal conditions. This point 
will be pursued in Section 3.5. 



3.4. An alternative approach to the conditions for dis- 
tortionless re-sampling 

The analysis presented so far has used the frequency 
domain approach. This is easy to use since the effects of 
cascaded linear processes may be compounded by simple 
multiplication of their frequency responses. However, the 
meaning of the results is not easy to visualise. In this 
Section we will examine the problem from the viewpoint of 
the time domain. 

Consider the system of Fig. 4{b) once more. The 
first sampler produces an output, s 2 (?), which may be 
thought of as a series of impulses whose amplitudes repre- 
sent the sampled values of the input signal s t (t), i.e. 



s, (t) = 



2 



a„S[t-nT ) 



where 8{t — t) denotes a unit impulse at t = t and a n = 
s,(«r ). (12) 

This is of course an example of pulse-amplitude modulation 
(p.a.m.). The train of impulses passes through the 

arbitrary filter before resampling. No distortion is intro- 
duced if the original set of impulses is reproduced again by 
the second sampler, i.e. we require 



s 4 (t) 



E 



a„8(t-nT -T), 



(13) 



which is s 2 (t) delayed by T 

and where clearly T=r + kT Q , 

k being some fixed integer. 



(14) 



Consider the very simple case when a = 1 and all 
other a n 's are zero, i.e. 

s 2 [t) = 5 (t), a unit impulse at the time origin. 

The output of the filter is therefore its impulse response, 
g(t), thus s 3 (t) = g{t). This signal is sampled at instants 
t = nT + t to produce an output which we require to be 
simply the original unit impulse delayed by T, i.e. 



The impulse response g(t) is constrained to contain a 
series of zeros at all sampling instants except the displaced 
time origin T. This is the condition for no distortion of a 
single impulse, which has merely undergone a delay T. 
Any real signal with an arbitrary sequence of a n 's will 
similarly have each of its impulses simply delayed and this 
is also undistorted provided the constraint on g(t) is 
satisfied. 

We may translate (15) into the frequency domain. 
Noting that the R.H.S. represents sampling g(t) at t = 
t + nT Q we may use the result of Section 7.1 for its 
spectrum. Thus, taking the Fourier transform of both 
sides we obtain 



-jo) T . 



z 



G{f-nf )e-l nb >° T 



(16) 



If we apply our assumption that Gtf) = for |/| >/ we 
obtain for the frequency range </</ 

e-> wr = G(fl+e" /w ° T G(f-/ ) 

This is identical to Equation (7) derived in Section 3.3. 
As before, the equation for negative/ is also easily obtained. 

Thus by argument in the time domain we have 
obtained the same condition for distortionless re-sampling 
as obtained earlier by a frequency domain approach. In so 
doing we have uncovered the meaning of this condition in 
the time domain, namely that the impulse response of the 
filter, g(t), shall be zero at all sampling points except the 
displaced time origin T. Each impulse resulting from the 
first sampling operation is blurred out in time by the filter, 
but because of these zeros \r\ g(t) there is no crosstal k f rom 
one sample to another at the points of sampling. Thus the 
original sample values are recovered by the second sampler. 

In fact our result is already well known as the con- 
dition for zero intersymbol interference in telegraph signal- 
ling, first described by Nyquist in 1928. 5 We may think 
of the series of impulses, s 2 (t), leaving the first sampler as 
being a series of telegraph symbols. The filter Gif) then 
represents the telegraph transmission line which delivers 
a distorted version of the series of symbols. However, by 
sampling at the appropriate instants the correct symbols may 
be decoded provided that there is no intersymbol inter- 
ference, which requires that the condition we have derived 
is satisfied. 



8{t-T)=g(t) 



E 



5{t-nT -T) (15) 



n = _oo 



Thus we have a constraint on the values of g(t) at instants 
t = T +nT a , namely g(T) = 1 

and g(T + mT ) = for al I values of m except m = 0. 



The analogy with telegraphy is still closer when we 
consider the actual sub-Nyquist coding system. Each 
sampler is then followed by a quantiser, the whole con- 
stituting an analogue-to-digital converter (a.d.c.) so that 
each sample is described by, say, / bits. A d.a.c. reproduces 
each sample which is thus constrained to one of 2 levels. 
This now corresponds to 2 -level telegraph signalling, con- 
veying/bits of information per transmission symbol. 

As a final illustration in the time domain we include 
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Fig. 7 - Illustration of the properties of the impulse response of a pair of comb filters 

{a) Assumed frequency response of a pair of comb filters ib) The resulting impulse response, showing sampling points 



two Figures which bring out the point that the impulse 
response of an antisymmetric filter of the form under dis- 
cussion is zero at all sampling points except the one at the 
(displaced) time origin. The Figures* were produced by 
means of a computer programme which performs Fourier 
transformation, the results being drawn by a plotter under 
computer control. 

Fig. 7 shows the frequency response and impulse 
response of a comb filter pair. The comb spacing chosen is 
coarser than for the filters actually used in order to give a 
clearer picture, but the important point is that the sampling 
frequency (396 arbitrary units) is still (n + 1 / 2 ) times the 
comb spacing (8 units). 198 units corresponds to / /2; 
note that the frequency response is antisymmetric about 
this point, at which it has the value V 2 . The sampling 
instants are superimposed on the impulse response and it 
may be seen that the value of the impulse response is zero 
at every sampling instant except the origin. This shows 
why resolution is not progressively lost as sub-Nyquist 
sampling systems using comb filters are cascaded. Each 
comb filter does average h.f. information from two suc- 
cessive scanning lines, but which the resulting signal is 
sampled this averaged information is discarded. 

Fig. 8 shows the frequency response and impulse 
response for a typical smooth antisymmetric low-pass res- 
ponse. In this case the impulse response extends over a 
shorter time span but retains the property that it is zero at 
every sampling instant except the origin. 

* The authors are indebted to A.H. Jones and R.E. Davies for these 
Figures. 



3.5. The effect of interposing a practical link 

In this Section we shall consider the effect of given 
imperfections in an intermediate analogue link when mixed 
digital and analogue transmission, of the form of Fig. 6(a) 
or (b), is in use. The main purpose is to see whether the 
performance specification of the analogue link needs to be 
any more stringent than that normally applied for analogue 
transmission. 

Start from Equation (5) which gives the 'effective 
insertion transfer function', Yif): 

Y(f) = Glf)+e->"° T G(f-f Q ) 

where G(f) is the combined response of the link and filters. 
Write G<f) = G x [f)G 2 if), where G l if) is the response of the 
filters and G 2 if) is the response of the link. Equation (5) 
thus becomes 

Yif) = G 1 (flG J (fl + e-i»° T G x \f- fo)G 2 (f- f ) (17) 

Suppose that the filters have an antisymmetric amplitude 
response A x if) with uniform group delay T x so that the 
condition for distortionless re-sampling (Section 3.3) would 
be satisfied if the link were perfect.* Specify that^ if) = 
for / > 5-5 MHz so that the link has to handle a signal of 
normal video bandwidth. This implies that A 1 if) is unity 
from to 3-37 MHz. 



* We have shown that this is true for either a pair of comb filters or 
the low-pass filter of Fig. 5. 
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(b) 

Fig. 8 - Illustration of the properties of the impulse response of a low-pass filter with frequency 

response antisymmetric about f /2 (= 198 units in this case) 
(a) Assumed low-pass frequency response (b) The resulting impulse response, showing sampling points 



For this range Equation (17) then becomes 

Y(f)=-\.G 2 (f) 

so that the insertion response in the band to 3-37 MHz 
equals the response of the link exactly in both amplitude 
and phase. 

We now pass to the h.f. response. To keep the 
analysis simple we shall assume that the link suffers gain- 
frequency variation or group delay variation but not 
both. 

Case (i) Link has uniform group delay T 2 but varying 
amplitude response A 2 (/). 

Substituting in (17) we obtain 



An interesting special case is when A l (f) has the 
smooth low-pass response of Fig. 5. As the cut-off about 
/ /2 (subcarrier frequency) is sharpened, the response of the 
link in the range above / Q /2 becomes less and less impor- 
tant. In the limit, A l If) becomes a perfect 4-43 MHz low- 
pass filter and the link performance above 4-43 MHz is 
quite unimportant. This result may appear somewhat 
startling until it is remembered that the ideal Nyquist 
signalling rate with 4-43 MHz bandwidth, i.e. / /2, is f 
symbols/second. The first sub-Nyquist sampling operation, 
sampling at / , has of course reduced the information 
content of the video signal to this rate. It can therefore be 
transmitted subsequently within a bandwidth / /2 without 
further degradation if ideal filters are used. 



Y(f) 



.4 1 (/M 2 (/) + { 1 - .4 , (fl J jyl 2 (f- / ) J 



-juiT 



! +7 V 



(18) 



provided that r, the sampling delay, is chosen to match the 
total group delay (T l + T 2 ). We are concerned with the 
sensitivity to departures of A 2 (f) from the ideal value of 
unity. If anything this sensitivity may be reduced com- 
pared with simple analogue transmission since the amplitude 
response \Y[f)\ at any frequency / between 3-37 MHz and 
5-5 MHz is a weighted average of the link response at fre- 
quencies / and if — f ), the weighting depending on the 
form of A x (f). For example, suppose a ±1 dB tolerance 
applies toA 2 if), The maximum variation in \Y{f)\ is still 
±1 dB but if the errors at /and (f— / ) are of opposite sign 
the variation in \Yif)\ will be reduced to less than ±1 dB. 



Case (ii) Link has uniform amplitude response but varying 
~ group delay 

We write G 2 (/) = e + - , *" ) where <pif) = phase response 
of link. 



Let 0(/) = <j>'lf) - coT c where T. 



-Mo/2) 



w /2 



thus0'(/ o /2) = O. 
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LY'(f) 



^(f-fj-tff)] 




reference axis $>(f)=-wT c 



-10 



Fig. 9 - Use of vector diagram to determine the overall response 

Substituting, Y{f)=A l (f)e~'" T c e^'^ + e-^° r (\-A x (f))e 

We note that the second sampler uses the subcarrier as 
a phase reference so that t = T c . The equation thus 
simplifies to 

Yif) = e->" T c iA 1 (/V ' (/) + (1 -A l (f))e>t'' {f - /o) j (19) 

The term within the brackets j |, say Yif), describes the 
deviation of Ylf) from unity gain and uniform group delay 
equal to T c . This is illustrated by the vector diagram of 
Fig. 9. 

The first observation to make is that the angle 
LY(f) (representing the departure of the phase response 
from uniform group delay) is less than or equal to the 
greater of (j)'(f) or <j>'(f — f ). This means that the phase 
response is no worse than that of the link. 

The amplitude response is represented by the length 
of the vector Y(f). At any given frequency this is reduced 



200 



when a link suffers from group-delay variation (Section 3.5) 
-/(w - o> )T C( j<l>'lf-f ) 



w 

c 
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5 5-5MHz 



3-37 4-43 



Fig. 10 - CCIR group delay template showing hypothetical 
worst responses as regards multiple coding 

( _'_ represent those group delay characteristics which give 

the greatest magnitude of 8(f)) 



from the ideal value of 1-0 as the angle \(p'(f- f ) - <j>'(f)\ 
increases. Thus the phase response of the link, 0(f), has 
given rise to a variation of the output in amplitude, as well 
as phase. To see whether this is within acceptable limits we 
need to know the likely variation of \<t>'(f — f ) —$'(/) J = 
0(/),say. 

Link specifications are frequently written in terms of 
group delay variation rather than phase response. Appen- 
dix 7.6 shows how 9(f) may be determined from areas under 
a group-delay curve. Fig. 10 is a group-delay template for 
a hypothetical reference circuit taken from the CCIR Green 
Book, Geneva, 1975, Vol. XII, p. 94. On it have been 
superimposed two curves which represent those character- 
istics giving rise to the greatest value of 6{f). (The response 
from to 3-37 MHz is immaterial for this purpose.) With 
either of these curves the calculated value of 9(f) for/ = 
3-37 or 5-5 MHz is 85°, using Equation (A11) of the 
Appendix. The minimum value of \Yif)\ is then given 
when A x (f) = A^(f - f Q ) = "A; \Y(f)\ is then 0-73, i.e. 
-2-65 dB relative to 1 -0. Suppose A t (f) takes the form of 
two comb filters in cascade. The minimum value of \Y(f)\ 
will then occur at frequencies where A l (f) = Vi, namely 
half way between peaks and troughs. These correspond to 
high frequency diagonal luminance information. The effect 
will only be significant towards the extremes of the combed 
range (3-37 and 5-5 MHz) because 9 -> as/-^/ /2. In 
contrast, if A x (f) is a smooth antisymmetric low-pass filter 
response the amplitude reduction we have described will be 
minimised since then 6 -+0 as A(f) -> %. 

Thus, if a practical link had a group-delay character- 
istic of the extreme form shown in Fig. 10 there would be a 
possible advantage in using the low-pass characteristic for 
A l (f)\r\ preference to the comb filter pair. 

We have considered amplitude and phase response 
variations separately; to consider both would require a 
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slight extension of the vector diagram technique illustrated 
in Fig. 9. 

In conclusion, we have shown that the introduction 
of an intermediate imperfect analogue link to a PAL sub- 
Nyquist system causes no drastic ill-effects. The insertion 
transfer function Y(f) resulting from such a link and the 
subsequent sub-Nyquist system obeys the following: 

(i) Y{f) is equal to G 2 (/) (the link response) at frequen- 
cies below 3-37 MHz. 

(ii) For frequencies above 3-37 MHz Y{f) depends on the 
response of the link at /and also at the 'mirror' fre- 
quency {f— f ), the two contributions being weighted 
according to the filter responses! ( if) chosen. 

(iii) For a link with varying amplitude response but uni- 
form group delay the maximum departure of Y[f) 
from unity will be no worse than the response of the 
link and may be rather better if A t (/) is carefully 
chosen. For example, if a link has poor h.f. response 
above subcarrier frequency then a low-pass filter with 
antisymmetric response cutting off fairly sharply 
above subcarrier frequency should be used. Y(f) will 
then show a better h.f. response than the link by itself. 

(iv) If a link has non-uniform group delay then Y(f) will 
show group delay variation and also amplitude res- 
ponse variation. When comb filters are used with 
such a link a small loss of diagonal information may 
occur but, for links satisfying the template of Fig. 10, 
this loss will be less than would be given by the comb 
filter without re-sampling. On the other hand, a low- 
pass filter will introduce a ripple in the h.f. response 
so that the information loss is smaller but more 
general. 

3.6. Multiple coding including quantising 

The discussion so far has considered repeated sampling 
operations in cascade. We shall now briefly consider the 
quantising that accompanies each sampling operation in a 
practical system of cascaded digital codecs. 

The assumption usually made is that when codecs are 
cascaded each contributes noise of mean-squared value 
equal to 5 /12, where 5 is the quantum spacing; 'noise- 
power addition' is thus assumed. As doubling 5 corres- 
ponds to using one less bit per sample we obtain the often- 
quoted result that using four codecs in cascade is equivalent 
to the loss of one bit per sample. 

This is not strictly valid. It rests on the assumption 
of a uniform probability density of the amplitude distri- 
bution of the signal which is quantised. This is perhaps 
reasonably true for the input to the first codec, but the 
signal applied to the succeeding stages may not necessarily 
have this property, especially if the sampling is in some way 
locked to the signal. 

Consider a special case in which: 



(i) the quantisers are identical and produce output levels 
V- t for inputs in the ranges V x — 5/2 to V\ +5/2 (where 
levels V\ are spaced by 5 ) ; 

(ii) any noise introduced between the quantisers is very 
much smaller than 5; 

(iii) sampling is in the same phase w.r.t. the signal in each 
case; 

(iv) the codecs are connected by a circuit of unity gain; 

and 

(v) an antisymmetric filter with uniform group delay is 
placed between the output of each quantiser and the 
input to the next so that Equation (11) of Section 
3.3 is satisfied. It is easy to see that the output of 
each sampler (excluding the first) will be the same as 
the output of the previous one. Thus no further 
quantising noise is introduced after the first sampling 
and quantising operation. 

In practice condition (i) could be relaxed somewhat 
to: each quantiser produces outputs V- t for inputs in the 
ranges V; - 5/2 to K, + 5/2 where V-, - 8/2<V-'< V } + 5/2, 
provided that the added noise (condition (ii)) is reduced 
indefinitely as the V t approach the extremes of the toler- 
ance specified. The quantising noise in this case is corre- 
lated; noise powers do not add and although the worst- 
case mean-squared noise is as much as 5/4 after only two 
codings, this remains the worst-case figure for an indefinite 
number of codings. This condition is more nearly prac- 
ticable: it would be satisfied by a perfect a.d.c. with a 
d.a.c. having a tolerance of ±5/2, or by an a.d.c. and d.a.c, 
each with ±5/4 tolerance. 

To analyse any further general relaxation of the con- 
ditions becomes very difficult, but it is possible to consider 
some other special cases, for example the subset of quan- 
tising noise known as contouring. This results when a 
slowly-varying signal is quantised; the output of the quan- 
tiser jumps from one quantum level to the next, spending 
a long period on each. 

Suppose that a sawtooth waveform is quantised. The 
result will be a staircase waveform with steps of height 5 
and having a regular spacing. If this signal is then re- 
quantised under the relaxed conditions given above the 
resulting contours will be irregularly spaced (in amplitude 
and position), but equal in number to those of the first 
quantising operation. If the gains are not accurately 
matched so that the quantum spacing of the second opera- 
tion, 5', is not equal to that of the first, 5, then a 'beating' 
effect takes place, illustrated in Fig. 11. Contouring is 
made more visible, but note that the figure shows an 
extreme case for clarity, in which the gain mismatch is 
large. 



4. Experimental work 

4.1. Introduction 

It was considered worthwhile to confirm the theore- 
tical results by a practical experiment in which a video 
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(a) 



ft) 



(c) 



Fig. 1 1 - The effects of multiple coding on contouring when there is a gain mismatch 

(a) Ramp waveform [b) The result after sampling and quantising {a) with quantum step S 

(c) The result after sampling and quantising (6) with quantum step 6 - "" '" 
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signal would be subjected to sub-Nyquist coding and de- 
coding several times in succession to see if picture quality 
would deteriorate further after the first operation The 
suggestion of using an antisymmetric low-pass filter as an 
alternative to the comb filter pair could also be tried. 

How was this to be done? A succession of sub- 
Nyquist codecs actually connected in cascade was imprac- 
ticable as this volume of equipment was not to hand. A 
system of recording and successive dubbing using only one 
codec was thus needed. However, an analogue quadruplex 
video tape recorder (VTR) could not be used because: 

(i) the defects of analogue tape recording would build 

up with successive dubs, masking the effects we 
wished to observe; 

and (ii) two analogue VTRs would be needed in the inter- 
mediate stages, one recording and the other playing 
back. Only one machine was available. 

The solution lay in the use of an experimental digital 
recorder designed and built in the BBC Research Depart- 
ment. 

4.2. The experimental arrangement incorporating the 
digital recorder 

The experimental digital recorder can record PAL 
video signals on magnetic tape in digital form. As organised 
at the time of the experiments it used the sub-Nyquist 2/ sc 
sampling system, with 7-bit p.cm. coding. Facilities are 
provided for error correction and concealment. It has 
separate record and play heads and associated circuits, so 



that in normal use play-after-record off-tape monitoring is 
provided after the fashion of normal studio audio recorders. 
This is a facility which is not provided by analogue VTR 
machines. 

When it is desired to examine the results of multiple 
processing the direction of tape motion is reversed so that 
the play head precedes the record head. Suppose the tape 
has an n\\ generation of processing recorded on it. This 
may be replayed, processed once more and recorded 'down- 
stream' of the play head as an {r + 1 )th generation record- 
ing. Thus an nth generation recording is produced by: 

(i) first recording from the picture source via the 
processing, 

and (ii) performing the dubbing and re-processing opera- 
tion (n — 1) times. 

The method is easily applied to our experiment since 
the same clock rate (2/ sc ) is used by the experimental pro- 
cessing and the recorder. Fig. 12 shows the arrangement. 
The filter shown as G(f) was in fact either: 

(i) a post-coding comb filter and a pre-coding comb 
filter in cascade, 

or (ii) a low-pass filter of approximately antisymmetric 
frequency response. 

The analogue VTR is not essential and was only used for 
convenience so that the results might be studied later with- 
out the digital recorder. Note that the final replay (stage 
three of the figure) uses no additional equipment to that 
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STAGE ONE! Initial recording 
pre-coding 
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input 



comb filter 




STAGE TWO! Dubbing with re-coding. (This stage is performed (/?-1) times) 
d.a.c. filter a.d.c. 
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STAGE THREE: Final replay 
d.a.c. 
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Fig. 12 - Method of producing an nth generation recording 



normally required for playback using the digital recorder, 
so the 'master' digital tapes may be viewed at any later time 
when the digital recorder is available. 

Use of the digital recorder for this experiment has the 
following advantages and disadvantages: 

(i) No analogue imperfections are introduced by the 
multiple recording process. The filters, a.d.c. and 
d.a.c. used are all part of the processing experiment. 

(ii) Only one machine and one codec are required for 
multiple recording. 



(iii) Any impairments introduced by the recording process 
take the form of digital errors which are easily dis- 
regarded since their visible appearance is quite dis- 
tinct from the imperfections we are studying (e.g. 
resolution loss). 



(iv) The digital recorder only has the capacity to deal with 
7-bit p.c.m. samples. Ideally we should study the 
cascading of 8-bit sub-Nyquist codecs; the restriction 
to 7 bits emphasises the effects of quantising, viz. 
contouring and quantisation noise. 
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Fig. 13 - A transversal filter having antisymmetric frequency response 
(a) Block diagram (b) Impulse response 



4.3. The filters used for G(f) 

(a) Two comb filters 

The pair of post-coding and pre-coding filters which 
were used were of the form described in Ref. 3 and were the 
ones used in stages one and three. Note that the total 
expression for G[f) with this combination must include the 
effect of the 5-5 MHz low-pass filters included in the a.d.c. 
and d.a.c. 



(b) Antisymmetric low-pass filter 

Appendix 7.5 demonstrates that any transversal filter 
of the form shown in Fig. 13 must have a frequency res- 
ponse which has the antisymmetric form required to satisfy 
Equation (1 1 ). This means that it is easy to design a filter 
having such a response using this configuration; an LCR 
filter would be more difficult to design. 

The overall filter transfer function G[f) referred to in 
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the analysis of Sections 3.2 to 3.4 must of course include 
the effect of the 5-5 MHz low-pass filters in the a.d.c. and 
d.a.c, also any filters which might be present in an analogue 
link. We must therefore choose the coefficients a r of the 
transversal filter so that its response above 5-5 MHz is 
negligible - its response will then dominate so that the 
overall response satisfies the antisymmetry condition within 
a sufficient accuracy. Note that ripples in the range 3-37 
to 5-5 MHz are entirely acceptable because they will auto- 
matically satisfy the antisymmetry condition. However, 
any ripples in the range to 3-37 'MHz will appear in the 
final insertion response since their 'mirror counterparts' lie 
above 5-5 MHz and are removed by the extra low-pass 
filters. 

It was decided to make a transversal filter to test the 
ideas in practice. To keep it very simple, a target specifi- 
cation appropriate to only one passage of the signal through 
it was chosen, namely that the ripple in the range to 
3-37 MHz was restricted to ±1 dB. Calculations showed 
that this could be achieved by a filter having a total of 5 
taps, i.e. 2 coefficient values a 1 and a 2 in addition to the 
pre-determined value of 1 / 2 for the centre-tap coefficient. 
In fact, a x and a 2 were allowed a range of adjustment so 
that some compensation could be made for the. short- 
comings of the delay elements used by departing slightly 
from the ideal settings. 

4.4. Experimental results 

Multiple recordings were made using several slides and 
one programme sequence. In every case the conditions of 
1, 2, 3, 4 and 5 sub-Nyquist coding operations were 
recorded. Two sets of recordings were made: one used 
the comb filter pair in the intermediate re-coding stages, 
the other used the transversal low-pass filter. 

No build-up of diagonal resolution loss or vertical 
colour resolution loss was observed, as we have predicted. 
However, some imperfections did build up as follows: 

(i) It was difficult to ensure that the gain of analogue 
parts of the system (between d.a.c. and a.d.c.) was 
exactly unity, so a slight change of gain and 'sit' was 
observed. 

(ii) A gradual build up of quantising noise and contouring 
occurred which was more noticeable on some slides 
than others. This was to be expected and was a result 
of being limited to 7-bit accuracy by the digital 
recorder. 

(iii) Using the comb filter pair in the intermediate stages 
a very slight loss at low frequencies occurred which 
caused visible streaking following vertical luminance 
edges. Section 3.5 shows that the response Y{f) at 
low frequencies follows G{f). From this we may 
deduce that the observed l.f. loss in the overall res- 
ponse was merely due to a slight droop in the comb 
filter response at l.f. which was magnified in effect by 
the multiple passage of the signal through it. Note 
that the last stage reached, i.e. 5 codings and de- 
codings, involved 5 passages of the signal through the 



comb filter pair so that effectively a total of 10 comb 
filters was involved. 

(iv) Using the transversal filter in the intermediate stages 
produced two defects. One was 'ringing' and em- 
phasis of vertical edges i.e. ripple in the frequency 
response. This was not surprising — the original 
filter was designed so that even in the ideal case it 
had ±1 dB ripple in the range to 3-37 MHz. The 
practical filter was slightly worse owing to imper- 
fections in the delay sections. Thus, when the signal 
passed through this filter several times (4 times for 5 
codings) this ripple was quite visible. 

The second defect was relative delay between chromi- 
nance and luminance i.e. difference in group delay for 
low and high frequencies. This is also 'caused by the 
delay elements which exhibited some dispersion. 

In summary, apart from quantising effects which are 
inevitable in any multiple digital coding scheme the only 
degradations introduced are caused by imperfections in 
analogue elements — filters, etc. — and were of the form 
which would equally be caused in analogue systems con- 
taining such elements. In particular, the sub-Nyquist 
nature of the digital coding used did not introduce any 
special restrictions on the use of multiple coding. 



5. Conclusions 

It is possible to cascade sampling systems without 
additive degradation due to sampling provided that certain 
conditions are met. In particular, the sub-Nyquist sampling 
of composite PAL video signals may be repeated indefi- 
nitely without any build-up of fundamental degradations 
introduced by the comb filters used in this technique. 

The possible future use of mixed analogue and digital 
transmission does not present a serious obstacle to the use 
of the sub-Nyquist technique for PAL video signals. An 
alternative filter for use at intermediate coding points is 
specified which is somewhat simpler than a comb filter and 
may permit a less stringent tolerance of the performance of 
intermediate analogue links. 
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7. APPENDIX 



7.1. Allowance for delayed sampling 

Consider a sampler which samples an input signal 
x(t) to produce a sampled output x s {t), Fig. 14(a). If the 
sampling occurs at instants t = + nT Q , classical analysis 
gives the well-known result for the output spectrum: 



X 5 (f) 



E 



X[f— nf Q ), normalised gain being assumed. 



What if the sampling were performed at instants 
t = r + nT ? In other words, sampling occurs j later than 
in the case first analysed. Conceptually, the same sample 
values could be obtained if the signal were made to occur r 
earlier, sampled at t = + nT and then delayed by r again. 

Thus if x'(t) represents the negatively delayed version 
ofx(t) we may write 

X'{f) = e + 1 WT XV). 

Ux'(t) is sampled, at instants t= +nT , to gwex' s (t) we 
obtain: 



\^ lx<f-nf )<J {u - n »o)r\ 



If this sampled signal x' s (t) is delayed by r we obtain 
the final output, x" s (t), which is equivalent to the signal 
obtained by sampling x[t) at times t = t + nT Q (see Figs. 
14(6) and (c)). Thus 



X" s V)=e-'" T 



\* \x{f-nf )e? { »- n "° )T \ 



m 



X(f-nf )e->""° T 



(A1) 



7.2. Analysis of repeated-sampling system 

We will use the result of 7.1 to continue the analysis 
started in Section 3.2. The system under consideration is 
shown in Fig. 4(b). Use Equation (A1) to describe the 
action of the second sampler: 



X'M 



E 



X'(f-nf ) 



S A V) = 



£ 



SM-nf^e-^or 



(A2) 
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Fig. 14 - Analysis of delayed sampling 
(b) Delayed sampler (c) Hypothetical version of (6) incorporating (a) 
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for |/1 >f 



Noting that S 3 {f) = G{f).S 2 {f) and remembering our 
assumptions regarding low-pass filtering we may substitute 
from (3) of 3.2 into (A2) to obtain: 

S 4 (f) = G{f){S l {f)+S l {r + fo)+S l {f-f )\ + 

+ e->"° T . G{f-f ){s l (f-f ) +S l (f) +S l (f- 2/ )j 

+ e + i"° T .G(f + f Q )\S l if + f Q ) +S 1 (f+ 2/ ) +S l (/}} 

(A3) 

We may simplify this further using our assumptions that 

5^ = 0' 

G(f) =0 

Considering /positive, in the range </</ Q : 

S 4 (f)=G{f]\S l {f)+S l {f-f Q )\ + 

+ e->"° T G(f-f ){S 1 tf-f )+S l ty 

= {&/) +e->"° T Gtf~ f )\ ■ {s^ + S^- f Q )\ 

(A4) 

and for /negative, 0>/>-/ o : 

S 4 {n = Gif)\S l (f)+S l {f+f )\ + 

+ e +/ "° T .G(f + / ){S 1 \f + / ) +S l (f)\ 

= \G[f)+e + '^ T .G{f + f )\{s i {f)+S 1 <f + f )} (A5) 

7.3. Derivation and interpretation of Y{f) 

The effective insertion transfer function, Y(f), was 
defined by Equation (4) as 



Y(f) = 



S 4 if) 



Equation (3) relates S 2 (f) to S i if) and Equations 
(A4), (A5) above relate 5 4 {f) to S % (/). We may therefore 
obtain Ylf) by dividing Equation (A4) or (A5) by Equation 
(3), yielding the results: 

for f o >f>0: Y{f) = G{f) + e- ,u <>r Gif-f ) (A6) 

andfor0>/>-/ o : Y{f) = G{f) + e +/u> ° T G{f +f Q ) (A7) 

We will now show that Ylf) behaves like a real base- 
band filter. Such a filter has a real impulse response y(t). 
\f y{t) is real it follows directly from the Fourier transform 
definition that 

Y(f) = Y*(—f) and vice versa. (Y* denotes the complex 

conjugate of Y). 



Consider the case of positive / and writing / = —v, so 
that v is negative. Then 

Y(-f) =Y(v) 

= G{v) + e + '"° T G(v +f ) from (A7) (valid since v is 
negative) 

Y*[-f) = Y*(v) 



coniu- 



= G*{v) +e~i° J ° T G*(v+f ) taking complex 
gate of above 

= G(-v) +e~J"° T G(-v -f Q ) since G(f) describes a 
real filter 

= G[f) +e~ 1 '"'° T G(f-f ) by substitution 

= Y{f) by (A6) (valid since/is positive). 

This is the desired result so that y(t) must be real and 
the combined effect of the arbitrary, but real, filter G{f) and 
the second sampling operation is therefore equivalent to a 
real filter of response Y(f). 



input 




> output 



Fig. 15 - Simple comb filter 

7.4. Response of simple comb filter 

Consider the simple filter of Fig. 15. By inspection, 
its impulse response, h x tf) say, is given by 

*,(') = ysjsw + su-rj} 

Applying the Fourier transform 

H 1 {f) = % (l + 1.e - '' wr ij 

-jo)T l / +jojT 1 -fuTj) 



= e 



+ e 



We will demonstrate that this is true for the present analysis. 
(EL-127) -17- 



i.e. 



-jwT 1 

2 (uT, 
H l {f) = e cos I L 



If two such filters are connected in cascade the overall 
response H 2 (f), say, is given by 

H 2 U)-H i 2 {f) = e-i» T i cos 2 \^y\ 



■ T i 1+ cos co7\ 
-/wTjl l 



H 2 V) = e 



Thus H 2 (f) is of 'raised cosine' form with constant 
group delay T l . Note that for use in PAL sub-Nyquist 
sampling T t would equal 283% subcarrier periods so that 
T 1 = 567% T Q . For convenience we shall write this as 
Tj = (m+V2)T .* The amplitude-frequency response A{f) 
is the modulus of H 2 (f) 



.e. 



A(f)- 



1 + cos cor, 



and A{f-f Q ) = 



|l + COS(C0- COqITjI 



|l H-COSICO^J -2m+ 17T)| 

2 



1 - cos co7\ 



so that 



(1 + cos cor. ) (1 - cos cor, ) 
A{f)+A{f-f )= L+ L_ 



1. 



Thus Equation (11) of Section 3.3 is satisfied so that 
a pair of simple comb filters, with the specified delay ele- 
ments, may be used for Glf) in the arrangement of Fig. 4(b) 
without any degradation. 

7.5. Response of transversal filter of Fig. 13 

By inspection of Fig. 13 we can write down the im- 
pulse response of the filter as 



g(t) = a n b(t)+a n _ , 8(t-T Q ) + 



+ «,6{f-/t- ir ) + %fi(r-n- W -T ot2 ) 

* Such a relation (with a different value of m) holds for some other 
sub-Nyquist systems. 



+ a^(t - nT ) + a 2 b(t - n + \T ) + 



+ a n _ ^8(t- 2n-2T )+a n 8(t- 2n - 1T ). 



le.g(t)=VS[t- 



8(t-n-rT n 



r= 1 



+ 8(t-n+r- ir o ) 
Taking the Fourier transforms of both sides we obtain 



G (/) = i/ 2 . e -/w< , /= + "- 1)7*0 + 
n 
+\^ a r (e~' 



ju>n - rT 



o + e -l'">n + r- 1T ) 



r= 1 



-ju 



2«-1 



♦E« 



%+ > a r \e 



r= 1 



2/--1 



+ e 



2/--1 

-/"I /'0 

2 



<2n--\ 



% + 2 } fl r cos ( )co:T D 



2 D 



/• = 1 



We see that the filter has uniform group delay (2n — 1 )/2T 
which is half the total of the delay sections used within it. 
The amplitude frequency response is 



A{f) = Vi + 2y a r cosl — )cor o . 



£' 



r = 1 



Consider the rth term, -4-1/} say, where 



A r (f) = a r cos | 



!2r- 1 



\o>T n 



thenA r {f) +A r (f-f ) = a r cos 



<2r- 1 



« r o + 



2 

/2/- - 1 
+ cos l J ( w _ oj )T \ 
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(2r-1)(2<0-w o )7' o (2r-1)co o r o 

7a r cos . cos 



= since cos (2r - 1 )- = 0. 
2 



Thus.4 (f)+A(f- f ) = V, + % + 2 j)T|4,{/) +A r {f- f )\ 



r = 1 



= 1. 



The filter therefore satisfies the condition for anti- 
symmetry of its amplitude-frequency response. 

7.6. Derivation of phase response 0(f) from group delay 
response Tif) 

Suppose a network has phase response 0(f) and group 
delay response T(f). They are related by 



Tif)-- 



-d<f>(f) _^d$ 
dco 2n df 



(A8) 



Integrating, we obtain 



£ 



0(f)- 0(0) = -2tt 



Tif) df 



i.e. 0(f) = — 2w 



Tif) df, since 0(0) = for a real network. 

(A9) 



4 ~' :i '<&. 



-f n 



/2 f-f 



Suppose Tif) is specified as a nominal group delay T with 
variation ATif), i.e. Tif] = T + ATif). 



f 



Then 0(f) = -2?r 



Tif) df 



-uT -2n 



ATif) df (A10) 



We may apply this result to the calculations of 
Section 3.5. In this case, we wished to know the value of 
6(f), where 



0(f) = 0'(f-f o )-0'(f) 



and 0'(f)= 0(f) + uT c where r c 



-0(f o /2) 
w /2 



Substituting, 

0(f) = 0(f- /„> + < w - w o> r c - 0V) - <^c 
= 0(f-/ o )-0tn-co o r c 
= 0(f-f o )-0(f)+20(f o /2) 
Now apply Equation (A10): 

f-fo 



6if) = -ioo-oj )T-2Tr 



ATlfldf 



+OJT + 27T 



ATif)df + 



[PTO] 



%± 



.frequency 



/2 f 



(EL-127) 
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foil 



-2.G) /2.T - 2n 



AT(f)df 



-2tt 



f-fn 



fo/2 

r 



AT(f)df+ 2 



■ 



AT{f)df- ATifldf 



Noting that ATlf) = AT{-f) for a network with a real 
impulse response this becomes 



m = -2n 



+ /o/2 / 

AT{f)df- ATlfldf 

f-fo 



-A 



0/2 



Refer to Fig. 16. The value of the term within the brackets 
[ ] is given by 

[ ] = (area A + area C) — (area B + area C) 



= area A — area B. 
Thus 0(f) = 2tt (area B - area A). 



(A11) 



SMW/SB 
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